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Abstract 

The decay K~ — > n^eu has been studied using in-flight decays detected with the 
"ISTRA+" spectrometer operating in the 25-GeV negative secondary beam of the U-70 
PS. About 550K events were collected for the analysis. The parameter of the vec- 
tor form-factor has been measured: A+ = 0.0286 it 0.0008(stat) ± 0.0006(syst). The 
limits on the possible tensor and scalar couplings have been obtained: /t//+(0) = 
Omito°o7t{stat)±0.026{syst); fs/f+{0) = 0.002l:[{:[J|^(stat) ± 0.003(s?/st) 



1 Introduction 

The decay K — > ez/7r°(Ke3) is known to be a very promising one to search for an admixture of 
scalar (S) or tensor (T) contributions to the Standard Model (SM) V-A amplitude. This decay 
has been extensively studied over recent years in different experiments with the charged and 
neutral kaons. Some published results indicate the anomalous S and T signals [1, 2]. 

On the other hand, a recent KEK experiment [3, 4] with a stopped beam and our 
preliminary studies [5] do not observe any visible contributions of tensor and scalar interactions. 
Another goal of our study is a precise measurement of the V-A f+{t) form- factor of the K^^i 
decay which is interesting in view of new calculations in the Chiral Perturbation Theory(ChPT) 
up to the order [6] . 

We present a new study of the K~.^ decay based on the statistics of about 550K events. This 
anaysis is an update of our preliminary results [5] . 

2 Experimental setup 

The experiment has been performed at the IHEP 70 GeV proton synchrotron U-70. The 
experimental setup "ISTRA-I-" (Fig.l) was described in some details in our paper [5]. 




Figure 1: Elevation view of the "ISTRA-F" detector. 



The setup is located in a negative unseparated secondary beam. The beam momentum is 
~ 25 GeV with Ap/p ~ 1.5%. The admixture of K~ in the beam is ~ 3%. The beam intensity 
is ~ 3 • 10^ per 1.9 sec. of the U-70 spill. The beam particles arc deflected by the beam magnet 
Ml and are measured by BPCi -r- BPC^ proportinal chambers with 1 mm wire spacing. The 
kaon identification is performed by -t- C2 threshold C-counters. 

The 9 meter long vacuumed decay volume is surrounded by 8 lead-glass rings LGi -j- LG^ 
which are used as the veto system for low energy photons. The photons radiated at large angles 
are detected by the lead-glass calorimeter SP2. 

The decay products are defiected by the spectrometer magnet M2 with a field integral of 
1 Tm. The track measurement is performed by 2-mm-stcp proportional chambers (PCi^PCs), 
1-cm-cell drift chambers {DCi -7- DC3), and by 2-cm-diameter drift tubes {DTi DT4). Wide 
aperture threshold Cerenkov counters {C3 and C4) are filled with helium and are not used in 
these measurements. 
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The photons are measured by the lead-glass calorimeter SPi which consists of 576 counters. 
The counter transverse size is 5.2 x 5.2 cm and the length is about 15 Xq. 

The scintillator-iron sampling hadron calorimeter HC is subdivided into 7 longitudinal sec- 
tions 7x7 cells each. The 11x11 cell scintillating hodoscope is used for the improvement of the 
time resolution of the tracking system. MuH is a 7x7 cell scintillating muon hodoscope. 

The trigger is provided by 5*1 ^ 5*5 scintillation counters, Cq C2 Cerenkov counters, and 
the analog sum of amplitudes from last dinodes of the SPi : 

T = S*! ■ ^2 ■ 53 ■ 5*4 ■ Co ■ Ci • C2 • ^5 • S(S'Pi), 

where 5*4 is the scintillator counter with a hole to suppress the beam halo, is the counter 
located downstream the setup at the beam focus, S(S'Pi) requires that the analog sum of 
amplitudes from the SPi be larger than ~700 MeV - a MIP signal. The last requirement serves 
to suppress the dominating K ^ ixv decay. 



3 Events selection 

During runs in Spring (runl) and Winter (run2) 2001, 363M and 332M events were logged 
on tapes. This statistics is supported by about 260M MC events generated with Geant3 [7] 
Monte Carlo program. The MC generation includes a realistic description of the setup with 
decay volume entrance windows, tracking chambers windows, chambers gas mixtures, sense 
wires and cathode structures, Cerenkov counters mirrors and gas, the shower generation in EM 
calorimeters, etc. 




Figure 2: The ratio of the energy of the Figure 3: The 77 mass spectrum for the 
associated ECAL cluster to the momentum events with the identified electron and two 
of the charged track. extra showers. 
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The data processing starts with the beam particle reconstruction in BPCi BPC4. Then 
secondary tracks are looked for in the decay tracking system and events with one good negatively 
charged track are selected. The decay vertex is reconstructed by means of the unconstrained 
vertex fit of the beam and decay tracks. 

A clustering procedure is used to find showers in the SPi calorimeter, and the two- 
dimensional pattern of the shower is fitted with the MC-generated patterns to reconstruct 
its energy and position. 

The matching of the charged track and a shower in SPi is done on the basis of the distance r 
between the track extrapolation to the calorimeter and the shower coordinates (r < 3cm). The 
electron identification is done using the ratio of the shower energy (E) to the momentum of the 
associated track (P). The E/p distribution is shown in Fig.2. The particles with 0.8 < E/p < 1.3 
are accepted as electrons. 

The events with one charged track identified as electron and two additional showers in the 
ECAL are selected for further processing. The 77-mass spectrum is shown in Fig. 3. The tt" 
peak is situated at M^o = 134.8 MeV with a resolution of 8.6 MeV. 

The selected events are required to pass 2C K cutt^ fit. To minimize effects of a beam 
associated background and systematics the cut on the difference between the fitted Pk value 
and the mean beam energy (25.2 GeV for the first run and 26.3 GeV for the second) is applied: 
\Pk — Pbeam\ < 1 GcV (sce Fig. 4). The missing energy Ei, = Ek — E^ — E^^o at this stage of 
the selection is shown in Fig. 5 for run2 data. 




Figure 4: The beam kaon momentum after Figure 5: The E^ energy compared with 
2C Ke3 fit. The points with errors are the MC. The points are data, the dark 

run2 data and the histogram is MC. histogram is the background contamination. 
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The peak at low E^, corresponds to the remaining K~ — > tt^tt'^ background. The cut 
Eu > 0.3 GeV is apphed to suppress this background. The requirement that the 5C X — > 7r~7r° 
fit should fail is applied for further background suppression. 

The surviving background is estimated to be less than 1.5%. 



4 Analysis 

Performing the procedure described in the previous section, 112K and 440K events are selected 
in runl and run2 data. The difference in the event output is explained by the higher SPi 
analog sum threshold in the second run. The distribution of the events over the Dalitz plot for 
the run2 is shown in Fig.6. The variables y — 2Eg/MK and z — 2E-,^/Mk-i where E^ and E^^ 
are the energies of the electron and 7r° in the kaon rest frame, are used. The analysis of the 
MC-generated background shows that the background events are concentrated at the peripheral 
part of the plot. 

The most general Lorentz- invariant form of the matrix element for the K~ 1~v-k^ decay 
is [8]: 

M^^lY^u{p;){l+^'')[2mKfsA{PK + P.)J+HPK-P.)af-^^^ (1) 
2 niK 

It consists of scalar, vector, and tensor terms. The /± form-factors arc the functions of t = 
[Pk — P^kY ■ In the Standard Model (SM), the W-boson exchange leads to the pure vector term. 
The scalar and/or tensor terms which are "induced" by EW radiative corrections are neghgibly 
small, i.e nonzero scalar or tensor form-factors would indicate the physics beyond the SM. 

The term in the vector part, proportional to /_, is reduced (using the Dirac equation) to 
the scalar form-factor. In the same way, the tensor term is reduced to a mixture of the scalar 
and vector form-factors. The redefined vector (V) and scalar (S) terms, and the corresponding 
Dalitz plot density in the kaon rest frame {p{E^, Ei)) are [9]: 

p{E^,Ei) ~ A-\V\'' + B-Re{y*S) + C-\S\^ (2) 
^ = /+ + {mi/mK)fT 

A = mKi2EiE,-mKAE^)-m^{E,-]AE^) 
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B = mimK{2E^ - AE^); E^ ^ niK - El - E^ 
C = mlAE^; AE^ = - E^; E™"" = ^ ^ 



.2 _^2^^2 

2'mK 

The terms proportional to mi and mf can be neglected in the case of the Kes decay. 

Then, assuming the linear dependence of the /+ on t: = /+(0)(1 -|- A+i/m^) and real 

constants fs, fr, we get: 



2\2 



p{E^,Ei) ~ mK{2EiE,-mKAE^)-{l + X+t/mi) 
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(3) 



The procedure of the extraction of the parameters A+, fs-, fx starts with the subdivision of 
the Dahtz plot region y = 0.12 ^ 0.92; z = 0.55 ^ 1.075 into 40 x 40 cells. 



The number of events in a cell 
extraction of, for example, A+ and 



fs 



/+(0) 



of the Dalitz plot, in the case of the simultaneous 
is fitted by the function: 



W 



MC, 



hi) 



fs 



./+(0), 



(4) 



Here Wi are MC-generated functions which are built up as follows: the MC events are generated 
with the constant density over the Dalitz plot and reconstructed with the same program as for 
the real events. Each event carries the weight w determined by the corresponding term in 
the expression 3, calculated using the MC-generated "true" values for y and z. The radiative 
corrections according to [10] are taken into account. Then, Wi are calculated by summing up 
the weights of the reconstructed events in the corresponding Dalitz plot cell. This procedure 
allows one to avoid the systematic errors due to the "migration" of the events over the Dalitz 
plot due to the finite experimental resolution. 
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Figure 6: Dalitz plot for the selected Figure 7: The ratio (A^^*^ - N^^)/Wmc 

K — * ez/7r° events. Run2 data. (see text) for the runl+run2 data. 

The parameters of the amplitude of the decay are extracted by minimizing the extended 
log-likelihood function: 
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where the sum runs over all populated cells of the Dalitz plot. Here "A" is the data/MC 

normalization parameter, Nf^^ - number of real data events in the (i,j)-cell, and N^^ - the 
MC estimated background. The background normalization is determined by the ratio of the real 
and generated K~ — > 7r~7r° events. The minimization is performed by means of the "MINUIT" 
program [11]. 

5 Results 

The results of the fit are summarized in Table 1. 

The combination of two runs is done by the simultaneous fit. The first line corresponds to 
the pure V-A SM fit. In the second line the tensor and, in the third, the scalar terms are added 
into the fit. 

The final errors are calculated by the "MINOS" procedure of the "MINUIT" program [11]. 
The errors are estimated at the 2 AC level of 2.3 that corresponds to the coverage probabihty 
of 68.27% for joint estimation of 2 parameters [13]. 





runl 


run2 


runl+run2 




0.0291 ±0.0018 


0.0285 ± 0.0009 


0.0286 ± 0.0008 


/t//+(0) 




0.024tH^^ 


0.02llH?^ 


/s//+(0) 


o.oost^i 


0.002lH^l 


o.oo2t°:°i 


xVndf 


1.13 


1.14 


1.14 


-^bins 


1120 


1119 


2239 



Table 1. Results of the fit 



To illustrate the quality of the fit. Fig. 7 shows the ratio of the real data minus MC-prcdicted 
background (A^^*^ — N^^) over Wmc versus t/m^o, calculated using expression 4, with A+ = 
0.0286 and /t = fs = 0. 

Different sources of systematics are investigated: 

• The Dalitz plot bining: 

AA+ = ±0.00006; A/t = ±0.003; Afs = ±0.0006. This kind of systematics is very small 
in our case because of the fitting method (see expression 4). 

• Selection cuts variation: 

1. The 5C X ^ 7r-7r° fit probability cut : 

AA+ = ±0.00041; A/t = ±0.0144; A/5 = ±0.0018. 

2. The 2C A' ^ ez/7r° fit probability cut: 

AA+ = ±0.00025; A/t = ±0.008; A/5 = ±0.0012. 
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3. The missing energy {E^,) cut: 

AA+ = ±0.00015; Afr = ±0.0081; A/5 = ±0.0008. 

4. The \Pk - Pbcaml cut: 

AA+ = ±0.00026; A/r = ±0.0104; A/5 - ±0.0018. 

5. The E/p electron selection cut: 

AA+ = ±0.00014; A/t = ±0.010; A/5 = ±0.0008. 

• The signal MC variation: AA+ = ±0.00005; A/r = ±0.001; A/5 = ±0.0001. 

• The background MC variation: AA+ = ±0.00015; A/t = ±0.011; A/5 = ±0.0003. 

From that, the total systematics is: AA+ = ±0.0006; Afr = ±0.026; A/5 = ±0.003. 

The comparison of our results with the most recent K"^ data [1, 4] shows very good agree- 
ment in the A+ parameter. We do not observe any visible contributions of scalar and tensor 
terms in the amplitude, in agreement with the conclusions of [4, 5]. 

some difference {2a) with the ChPT 0(j9^) calculat ions for A+: A+ = 0.031 [12], is observed. 

In addition, a possible contribution of the quadratic term A^t^/m^ into the /+ form-factor 
is searched for. The value of A+ = -0.00042t|]:|jg}^ is obtained with the fixed A+ = 0.0286. If 
A+ is allowed to vary, its value is shifted to A+ = 0.02867 ± 0.0020 and the quadratic coefficient 
becomes A^ = -0.002lHoi6- 

We conclude that a possible quadratic contribution into the 
vector form-factor is compatible with zero in both cases. 

6 Summary and conclusions 

The K^^ decay has been studied using in-flight decays of 25 GeV , detected by the "ISTRA+" 
magnetic spectrometer. Due to the high statistics, adequate resolution of the detector, and good 
sensitivity over all the Dalitz plot space, the errors in the fitted parameters are significantly 
reduced as compared with the previous measurements. 

The A+ parameter of the vector form-factor is measured to be: 



A+ = 0.0286 ± 0.0008 (stat) ± 0.0006 (syst). 




The limits on the possible tensor and scalar couplings are derived: 



/t//+(0) = 0.021: 
/5//+(0) = 0.002 



+0.064 
-0.075 



(stat) ±0.026 {syst); 
(stat) ± 0.003 (syst) 



,+0.020 
-0.022 
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